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This  thesis  expands  on  the  work  of  Myers  in  the  control  of 
linear  periodic  systeas,  it  applies  the  control  developed  for  the 
linear  system  using  aodal  variables  to  the  non-linear  equationfi  of 
motion.  While  only  one  example  is  investigated  in  this  thesis, 
the  technique  can  be  applied  to  periodic  systems  in  general.  This 
thesis  used  uncontrolled  modal  variables  with  scalar  control  and 
work  needs  to  be  expanded  involving  the  controlled  modal  variables 
and  other  control  schemes  to  fully  understand  the  effect  of  large 
perturbations  on  the  non-linear  equations. 
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AtWtMCt 


The  attitude  of  a  spinning  syanetrical  satellite  in  an 
elliptical  orbit  Is  analysed.  The  perturbed  notion  of  the 
satellite  is  described  by  linear  equations  with  periodic 
coefficients.  Stability  is  deternlned  by  Sloqnet  theory.  Active 
control  is  added  to  the  systen  and  results  lead  to  a  linear 
periodic  control  lav.  Scalar  control  fron  the  linearised  systen 
is  iaplenented  to  evaluate  the  perfotnance  of  the  control  law  on 
the  non-linear  equations  of  notion.  For  snail  disturbances,  it  is 
illustrated  the  controlled  non-linear  response  duplicates  the 
linear  case.  For  larger  disturbances,  phase  portraits  show  the 
resulting  behavior  of  the  coupled  nodes.  As  the  perturbed  notion 
increases  in  nagnitude  stability  regions  appear.  Initial  notion 
results  either  a  return  to  the  initial  equilibriun,  an  oscillation 
around  the  equilibriun  point  or  divergence  fron  the  initial 
equilibriun. 


is 


Chapter  1 
Introduction 


Since  satellites  vere  first  put  into  orbit,  designers  have 
been  concerned  with  Mintaining  the  attitude  with  respect  to  soae 
fixed  reference.  As  satellites  have  grown  Bore  complex  and 
expensive,  this  concern  has  increased,  tfith  this  Increased 
coaplexity  has  come  an  increase  in  the  need  for  precision 
pointing.  Attitude  control  of  a  satellite  is  now  of  extreme 
concern.  Through  the  years  various  methods  have  been  used  to 
achieve  this  attitude  control  (ref  1).  Satellites  in  near 
circular  oxbits  can  be  designed  to  take  advantage  of  the  gravity 
gradient  to  produce  attitude  control.  However,  this  method  has 
one  main  restriction.  Low  natural  frequencies  result  in  long 
decay  times  for  any  disturbance.  While  this  may  be  sufficient  for 
some  systems,  it  certainly  would  not  be  sufficient  for  any 
precision  surveillance.  Active  control  must  be  used  for  many 
satellites.  While  there  are  many  possible  control  systems,  each 
with  its  advantages  and  disadvantages,  the  one  of  interest  here  is 
the  gyroscope. 

The  gyroscope,  with  its  large  angular  momentum  and  torque 
free  motion  has  long  been  utilized  as  an  ideal  attitude  reference 
device.  The  next  natural  step  was  to  make  the  satellite  a 
gyroscope  by  spinning  it  about  one  of  its  axes.  The  stability  of 
a  spinning  unsymmetrical  satellite  in  a  circular  orbit  was 
Investigated  by  Kane  and  Shippy  (ref  5)  in  1963.  Kane  and  Barba 
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(ref  4)  also  studied  a  spinning  syaaetrical  satellite  In  an 
elliptical  orbit  in  1966.  Both  these  investigations  yielded 
equations  which  were  non-linear  and  non-autonoaous .  In  both  cases 
the  satellites  were  nodeled  as  a  tine  periodic  linear  systea,  with 
the  stability  deteralned  by  Floquet  theory.  In  Kane  and  Barba's 
study  of  a  satellite  in  an  elliptical  orbit,  it  was  deaonstrated 
that  the  stability  was  dependent  on  the  inertial  properties,  orbit 
eccentricity,  and  satellite  spin  rate. 

The  need  for  aore  precise  attitude  control  has  elialnated 
gravity  gradient  stabilization  on  these  satellites.  For  spinning 
satellites  in  an  elliptical  orbit,  many  configurations  are 
unstable.  Linearized  systems  with  periodic  coefficients  have 
recently  been  given  their  deserved  attention.  Hyers  (ref  7) 
investigated  this  type  of  systea  for  a  satellite  with  two  unstable 
Bodes.  The  scheme  employed  transforming  the  state  variables  to 
modal  variables.  Calico  and  Yeakel  (ref  2)  used  a  similar 
technique  with  one  unstable  mode. 

This  thesis  extends  on  the  above  work  by  taking  the  solutions 
of  the  linearized  equations  and  applying  these  solutions  to  the 
non-linear  equations.  Theory  states  that  for  'small '  disturbances 
the  non-linear  systea  should  be  consistent  with  the  linearized 
systea.  The  unresolved  question  which  remains;  How  small  is 
'small'?  The  modal  variables  plotted  on  phase  portraits  provide 
information  that  aid  in  determining  the  answer. 
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Chapter  2 
THBORY 


In  this  chapter  the  eqaatlons  of  Motion  of  an  spinning 
sywietrlc  satellite  In  an  elliptical  orbit  around  the  earth  are 
developed.  The  equations  describing  both  the  attitude  Motion  and 
the  trajectory  Motion  of  the  satellite  will  be  presented.  The 
trajectory  equations  are  uncoupled  froM  the  attitude  equations  but 
the  attitude  Motion  is  affected  by  the  periodic  orbital  Motion. 
The  attitude  equations  are  linearized  into  a  systeM  of  tlMe 
periodic  linear  equations  and  Ploquet  theory  Is  used  to  deteralne 
the  stability  of  the  systea.  In  Reference  1,  it  has  been  shown 
the  unstable  systea  aay  be  stabilized  by  using  feedback  control. 
This  control  aay  be  supplied  by  either  scalar  or  a  vector  control. 
Scalar  control  Is  developed  In  this  thesis.  Once  the  control  for 
the  linear  systea  has  been  developed  for  ’saall'  disturbances  the 
control  will  be  applied  to  stabilized  the  non-linear  systea. 


2.1  Equations  of  Notion 

The  development  of  the  equations  of  Motion  presented  below  is 
described  by  Kane  and  Barba  (ref  4).  This  thesis  develops  two 
basic  equations  for  the  trajectory  of  a  satellite  around  a 
spherical  syaaetrical  attracting  body.  As  described,  these 
equations  are  developed  in  any  introductory  astrodynaaics  text. 


i'  -  , 


r*ve  a*n  i^-e* 


(2.2.1) 


(2.2.2) 
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where; 


T  Is  defined  as  the  period  of  the  satellite.  The  orbital  eleaents 
r,v*,a,  and  e  are  defined  In  Figure  1, 


Figure  2.1  Slllptlcal  Orbit  Elements 


Equations  (2.1.1)  and  (2.1.2)  Is  non-dlnenslonallsed  by 
defining  the  two  non-dimensional  variables  for  time  and  distance 
as; 


T  «  nt  t»~ 
n 


C  -  I  r-ca 


(2.1.4) 


Solving  equation  (2.1.2)  for  u  yields; 

«  -  a*njl-e* 


(2.1.5) 
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Inserting  this  valae  in  equation  (2.1.1)  yields; 


r  4  n*a*  ^  - 

r* 


(2.1.6) 


Before  this  equation  can  be  non-diaensionalized,  differentiation 
with  respect  to  r  wist  be  defined,  (represented  by  priwss) 


"  d  fdr”)  d  fdTfdr'l'i 

*  dtldtj  arlit l3tJJ  •  ^  * 


Rearranging  equation  (2.1.6)  yields 


r  ,  («»-ll«V  ,  »V 

n*a  r*  a  n*  n*a  r* 


and  using  equations  (2.1.7)  reduces  this  to 

*  1  1 

C"+  5_Ii  ♦  i-  •  f 


(2.1.7) 


(2.1.8) 


(2.1.9) 


The  value  of  t^  or  can  be  arbitrarily  selected.  Therefore,  for 
convenience,  it  is  assueed  the  satellite  is  at  perigee  for  tBr^f. 
Hence, 


r  *  a(l-e) 

P 


C(B)  «  ^  -  1-e 


r  >  0 

p 


CMf)  *  f 


(2.1.18) 


Mon-dlnensionalizing  the  true  anonaiy  in  equation  (2.1.5),  and 
noting  prises  represent  differentiation  with  respect  to  t  produces 


du  6v  dt 

ar  '  fft  37 


1 

n 


(2.1.11) 


J 
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s 


d  (M  ,  d 

(2.1.12) 

. .-ijl-e*  ,, 

C* 

Mote  that  C  is  a  periodic  function  of  r  with  period  2n. 

To  develop  the  equations  of  Mtion  for  a  spinning  satellite 
in  an  elliptical  orbit,  it  is  assosMd  the  satellite's  orbit  is  not 
effected  by  the  orientation  of  the  satellite  (i.e.,  the  satellite 
is  considered  to  be  a  point  mss). 

The  satellite's  orientation  with  respect  to  the  orbital  axes 
is  determined  by  defining  an  orbital  reference  frame  A.  Frame  A 
is  defined  such  that  A^  points  outifard  along  the  radius  vector. 
A^  will  be  perpendicular  to  A^  in  the  orbital  plane  and  defined 
such  that  at  perigee  and  apogee  the  satellite  will  be  moving  in 
the  positive  A^  direction  only.  A^  will  be  perpendicular  to  both 
A^  and  A^  in  such  a  way  to  form  a  right  handed  coordinate  system 
(illustrated  in  Figure  2.1). 

The  body  axes  is  defined  in  frame  X  and  is  obtained  from 
frame  A  by  a  l-2>3  rotation  through  angle  0,0  .  and  B 

ft  *  9 

(illustrated  in  Figure  2.2).  The  satellite  is  symmetrical,  having 
two  equal  moments  of  inertia,  therefore,  a  nodal  axes  will  be  used 
which  is  one  rotation  away  from  the  body  axes.  Thus  it  can  be 


C* 
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obtained  through  a  1-2  rotation. 


Figure  2.2  1-2-3  Rotation 


The  angular  velocity  of  the  orbital  reference  fraae  k,  with 
respect  to  an  inertial  fraae,  and  the  angular  velocity  of  the 


satellite,  expressed  In  the  nodal  axes  are  respectively. 


-A/'l  • 

<0  s  ve 


-X/a 
0)  ■ 


(l> 


Ci> 


(•> 


(2.1.13) 


Stepping  through  the  rotations  will  develop  the  orbital  rates  of 
the  different  fraaes,  . 

Rotations  - 


1)  about  through  9^  to  (Fraae  B) 


-■/A  i  ^  A  ^ 

“  •  V./ 


11)  about  B.  through  9  to  F^  (nodal) 

m  m  G 


-OTi  A  <">  X  ^ 

“  ■  V./  ®.'c. 


(2.1.14) 


(2.1.15) 


The  angular  velocities  aay  be  added  %dien  In  the  sane  reference 
frane. 
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n  •  n  — 

The  transforaatlon  aatrices  used  in  equation  (2.1.16)  are  defined 
an  follovsf 


■  cos©^  t  - 

■8ln©^‘ 

L  s 

CB 

f  1 

t 

(2.1.17) 

sin©^  • 

m 

cos©^ 

'  1  0 

f  ' 

L  * 

aA 

f  COS©^ 

sln©^ 

(2.1.18) 

1  -sin©^ 

cos©^ 

cos©^ 

ain©^ain©^ 

-cos©^8in©^' 

• 

cos©^ 

sin©^ 

(2.1.19) 

sin©^ 

-8ln©^C08©^ 

cos©  cos© 

<  s 

Bxpandinq  equation  (2.1.16)  and  uninq  equations  (2.1.17)  and 
(2.1.19),  yields  the  expression  tor  angular  velocity  expressed  in 
frane  C,  the  nodal  fraae. 


cxxl 
'  )  * 


f 

1 

f 


-sin© 

cos© 
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s©^  sin©^sln©^  -co8©^8in©  1  f  •  1 
•  cos©^  sin©^  1  1  *  I 

n©^  -8ln©^C08©^  C08©^C08©^|  I  V  I 


.1.29) 


Perforaing  the  linear  algebra  yields; 


f  C/a 

1“ 


}■ 


♦  -V  co8©^8in©.?  4  u  8ln©^e_  4  v  co8©  cos©  e 
t  a  1  t  s 


(2.1.21) 


Or  written  in  aatrix  fora; 
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0)  i  • 

Jc 


0^coa6^~  u  cosd^slnd^ 
oJ+  V  sln®^ 

3  8lii©,+  V  coad  cosd 

S  mm 


(2.1.22) 


This  Is  the  angular  velocity  vlth  respect  to  the  Inertial  £raae 
written  In  the  nodal  fraae. 

Ill)  The  last  rotation  0^,  la  about  to  obtain  to  the  body 
fraae  X.  It  la  expressed  in  the  nodal  fraae  In  order  to  achieve 
the  ilnal  fora  as 


r  M  ^ 


{• 


f  y 

e. 


(2.1.23) 


(2.1.24) 


n 


9^eo9e^-u  cosG^sind^ 
0^4  V  sinO^ 

d^alnd^+v  co8®^co8d^+0^ 


(2.1.25) 


Defining  «^s  9^  ox, 


produces. 


<ji'  •  9  8in0  *  V  coa9  coa9 

fl  A  S  AS 


(ij 


(2.1.26) 


(2.1.27) 


Angular  aoaentua  vector  about  the  center  of  aass  Is  defined  as; 

{ "ci/OJc  n 

Assualng  a  syaaetrlc  satellite,  the  aoaent  of  Inertia  aatrix  has 
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the  £ora; 


therefore; 


{"=1 


■  A 

0 

0  ‘ 

’  u 

0 

A 

0 

<•> 

s 

0 

0 

C 

<l> 

•J 

From  Euler's  Moaent  equation; 


^  "=}  ■  { "=1 


(2.1.29) 


(2.1.31) 


(2.1.31) 


«fhere 


M  "4  ■  { "4. 


'“Cy'I 

(j  ■ 


•  « 

«;  I 


(2.1.32) 


(2.1.33) 


Noting 
becoaes; 


In  equation  (2.1.28)  is  constant,  equation  (2.1.32) 


A  0  f 

0  A  f 
0  0  C 


<*> 


0  ^ifA  0  0if«'|  f  Ml 

«;  0  -«  0  A  0  j  «  H  ! 

"i  •  J[*  •  ^  Jl"sJ  I**.] 


(2.1. 


Again,  linear  algebra  operations  yield; 


A«, 

C«. 


♦ 


0  «  1  f  A«  n  f 

«;  0  -«  I  ]*«,[-  j 


(2.1.35) 


Further  linear  algebra  gives  three  equations; 


34) 
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(2.1.36a) 


Jtoj*  C<o^«,  ■  (2.1.36b) 

(3»^-  kbi^ta^-k  «  M^  (2.1.36c) 

Expanding  equations  (2.1.36)  with  the  definition  of  gives 


Ka  -JUd  <d  4  AA)  ©  -f  Cb>  «d 

(2.1.37a) 

i  so  s  a  s  a 

A 

•  • 

Au  -JUt)  (i>  -  Acd  ©  -  Cfa>  (d 

*  K 

(2.1.37b) 

a  i  a  A  a  t  a 

a 

ci. 

•“a 

(2.1.37c) 

Dividing  by  the  Moaent  of  inertia  k,  and  coabining  terns 

produces 

•  C-A  *** 

<d  ♦  <d  <d  4  <d  ©  ■  — 

«  A  as  as  A 

(2.1.38a) 

•  A-C  \ 

61  ♦  -  ■'  6)  6>  ^  61  d  » 

«  A  t  ft  Aft  A 

(2.1.38b) 

M 

a  C 

(2.1.38c) 

Vith  the  applied  aonents  being  (4;4f3) 

1 

II 

a'* 

(2.1.39a) 

Mjj«  3ilj  (A-C)  sln©^cos©^ 

(2.1.39b) 

V* 

(2.1.39c) 

To  sinplify  notation,  define; 

Km^ 

A 

(2.1.48) 

and  rewrite  equations  (2.1.38)  as; 

<d^4  K<d^<d^4  “a^a*  * 

(2.1.41a) 

-3n*Kain®  cos© 

- 5 -  (2.1.41b) 


•  (2.1.41c) 
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Taking  the  derivative  with  respect  to  tlM  of  equation  (2.1.25) 
also  results  In  expressions  for  the  rate  of  change  for  the  angular 


velocity, 

u  »  S  cos6  -  i^osd  slnd  ♦  vO  slnd  sin©  -v©  cos©  cos© 

S  AS  •  ft  SS  ft  z 

©^*  vsln©^-l>  i^^cos©^  (2.1.42) 

«  «  ©  sin©  ♦  ©  ©,sln©,+  vcos©^cos© -v©^sln©  cos©  -w©  cos©  sin©  +© 

ZtL  ZAZ  Z  AZft  ft  ZZ  ft  zz 

The  next  step  In  solving  the  equations  of  notion  Is  to  linearize 
the  equations  about  an  equlllbrlun  point.  Later,  this  solution 
will  be  Inserted  into  the  non-linear  equations.  The  equillbrlun 
point  Is  chosen  such  that  ©^«  ©^«  f.  ©^  Is  chosen  after  snail 
angle  approximations  are  made  on  equation  (2.1.22); 


(2.1.43a) 

MB©  ^  vB 
s  s  t 

(2.1.43b) 

»  ■  ©  *  V 

m  a 

Integrating  equation  (2.1.43c)  yields. 

(2.1.43c) 

«t»©  +  ©  *  U  f  V 

s  a  s  o 

o 

(2.1.44) 

dropping  the  subscript  on  <«>  and  with  ©^  ^  f, 

o 

solving  for  ©^; 

©^«  Mt-v 

(2.1.45) 

where  ««  ■  n(l4«()  with  n  defined  In  equation  (2.1.3)  and  oi  Is  the 
constant  spin  rate  about  the  axis.  With  these  values 
substituted  in  equations  (2.1.42)  and  (2.1.39),  the  accelerations 
f  and  the  external  moments  t.  With  this  equilibrium 

point,  the  axis  remains  perpendicular  to  the  orbital  plane.  If 
this  Is  a  stable  equilibrium  point  the  satellite  will  return  to 
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this  orientation  if  disturbed  by  a  'smII*  aaount.  To  deteraine 
if  this  is  a  stable  point,  state  variables  are  defined; 

K  K  ]  (2.1.46) 

and  the  perturbation  of  the  state  vector  in  the  linearized  systea, 

x'»  te  (2.1.47) 

The  first  two  rows  of  the  A  aatrix  are  trivial,  the  reaalnlng  two 

rows  are  found  by  taking  derivatives  of  equations  (2.1.43); 

<|>^«  (2.1.48a) 

<1)^»  e^+  ve  (2.1.48b) 


«  •«  •• 


(2.1.48c) 

Substituting  the  values  of  <a  froa  equations  (2.1.43)  and  w  froa 
equations  (2.1.48)  in  equation  (2.1.41a)  yields; 

J  *  (2.1.49) 

*  •'^eJV  ■  K"  (2.1.58) 

and  using  the  chain  rule,  the  derivative  w.r.t.  r  is  produced, 

n  •  -nv  ln®^4  Knv  ♦  "  **  d^+lnv  -  ne^-  Kn|»>  11"®, 

This  reduces  to  the  identical  equation  as  (2.1.51)  but,  w.r.t.  r; 

*♦  ***  ♦  '*  [**  ■  ®s"  *^1**  ♦®,JJ®,<2.1.51) 

Equation  (2.1.51)  produces  the  third  row  of  the  A  aatrix. 
Utilizing  equation  (2.1.41b)  and  again  using  saall  angle 
approxiaations  on  the  right  side  produces  an  equation  for  d^; 
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-3n*Ke 


(v^O  (V“)  ■  '*•*•”’ 

Changing  the  derivatives  to  primes,  solving  for  and  factoring 


teras  yield  the  fourth  row; 


.1.53) 


Bquation  (2.1.47)  may  now  be  trritten  as, 
d'  1  [  I  «  1  f 

A 

b'  I  I  •  1 

s 

.  ■ 

®r  *  *. 

®i'  "  *«  * 


%fhere  elements  of  A  are  defined  as  shown  and; 


(2.1.54) 


-V' 

V'  ^  ®e) 

.  ,  3H 

*..*  V  ^ 

A  is  now  expressed  in  teras  of  v' ,  v" ,  B^  and  C*  It  is  desirable 
to  have  A  expressed  in  teras  of  one  variable  C  and  its  derivative 
C^.  This  can  be  easily  done  by  replacing  v'  and  v"  using 
equations  (2.1.11)  and  (2.1.12).  nay  be  eliminated  by  noting 
9,  hence  is  a  constant,  or  from  equation  (2.1.43c), 


v'*B_’‘  constant 


Let  the  constant  be  (o-fl),  then 


B‘  (a+1)  -  4i-4- 

s 


(2.1.55) 


(2.1.56) 


A  is  now  expressed  in  a  more  convenient  form  with  the  non-zero 


2-12 


eleaents  of  equation  (2.1.54)  as 


!• 


S4 


1 


I 


ss 


\ 


^.4“ 


— *  -  (a+l)(l+IC)Ep 

C  C 

-2fl  e*"  . 

C* 


- (a+l)(UK) 

C* 


4* 


3K 

C* 


2.2  PLOQUET  THEORY 

The  systea  here  is  tlae  periodic  with  period  T, 

A(t)«A(t+T)  (2.2.1) 

therefore  the  stability  can  be  found  using  Ploquet  Theory. 
Ploquet  discovered  this  technique  In  the  latter  halt  of  the 

llfi's.  Our  tiae  periodic  systea  is  represented  by; 

X'»A(t)x  (2.2.2) 

where  x  is  the  state  vector  and  X'  is  the  derivative  with  respect 
to  T.  The  variable  t  is  used,  instead  of  tiae,  as  a  non- 
diaensional  variable.  A(t)  Is  known  and  is  a  periodic  aatrlx  with 
period  T.  Any  solution  x(t)  aay  be  written  in  the  fora 

X(T)«y(T,T  )c  (2.2.3) 

O 

where  wlr.r  )  is  a  fundaaental  aatrix  and  c  a  constant  vector. 

O 
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Solving  for  c  yields, 


c»vr‘(T  ,r  )x(T  )  (2.2.4) 

o  o  o 

When  y(T  ,T  )  «  I,  the  identity  Matrix,  V'(t,t  )  is  teraed  the 

o  o  o 

Principal  Pundaoental  Matrix  and  (2.2.3)  becones 

x(t)»1(t,t  )x(t  )  (2.2.5) 

o  o 


I 


where  the  Principal  Fundaaental  Matrix  is  denoted  by  C  and  has  the 
property 


•(T  ,T  )■  1 
o  o 

►  By  letting  t^=8, 

•(S,i)«  I 


(2.2.6) 


(2.2.7) 


The  B(t,0)  matrix  has  columns  of  which  the  first,  is 

» 

obtained  by  solving 

^^»A(t)^^  (2.2.8) 


)  In  general. 


► 


or. 


(2.2.9) 


•  (T,|)«A(T)B(T,i)  (2.2.18) 
Equation  2.2.18  can  be  solved  numerically  using  the  initial 
conditions  from  equation  2.2.7.  Numerical  Integration  would  also 
be  used  to  simultaneously  solve  for  the  variables  C  ond  C' 
contained  in  the  matrix  A. 

All  of  the  above  equations  are  true  for  all  time  dependent 
matrices.  Ploquet  Theory  states  •(▼,8)  can  be  factored  into  two 


c 


€ 


€ 


Matrices,  P  and  J,  in  the  Manner, 

t(T,i)-P(T>e'*^p"‘(f)  (2.2.11) 

The  J  Matrix  is  a  constant  Matrix  and,  for  convenience,  is 
expressed  in  Jordan  noraal  fora.  The  diagonals  of  the  J  Matrix 
are  the  analog  of  the  eigenvalues  for  a  constant  coefficient 
systeM  and  are  called  the  PoincarM  exponents,  co^. 

The  Matrix  P(t)  is  a  tiae  periodic  Matrix  vith  the  saae 
period,  T,  as  the  original  systea  (equation  (2.2.2)).  If  the 
systea  had  constant  coefficients,  P  would  be  the  constant  Matrix 
of  eigenvalues.  Thus,  solving  the  Ploguet  problea  reduces  to 
deteraining  the  constant  Matrix  J  and  the  periodic  Matrix  P  over 
one  period. 

Evaluating  f (T,0)  (recall  T  is  the  systea  period),  yields 

1  (T,i)«P(T)e'*V*(i)  (2.2.12) 

However,  F(T)»P(I)  since  P(t)  is  periodic.  Therefore  equation 
(2.2.12)  becoaes 

#(T,i)«P(i)e''V*(i)  (2.2.13) 

The  Matrix  *(T,f )  is  called  the  aonodroay  Matrix  and  P(H)  is  the 

Matrix  of  eigenvectors  of  •(T,f).  The  technique  of  constructing 
the  aonodroay  Matrix  is  to  nuaerically  integrate  equation  2.2.11. 

Proa  equation  (2.2.13)  it  is  seen  *(T,B)  and  e*^^  are  siailar 

and  therefore  have  the  saae  eigenvalues,  \^.  If  the  were  real, 

J  would  be  a  diagonal  Matrix  with  the  Poincare  exponents,  on 

the  diagonal.  The  eigenvalues  of  J,  the  Poincare  exponents  o>^, 

JT 

are  related  to  the  eigenvalues  of  e  (or  ((T,!))  by 

Xj^sexpiWj^T)  (2.2.14) 


€ 
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or 

V  f  \  (2.2.15) 

If  the  K  are  coaplex  then  taking  the  log  requires, 

i  Jin  |XJ  +  i  arg(X^)j  (2.2.16) 

The  stability  of  the  system  is  now  available.  Since  P(t)  Is 
periodic,  the  stability  of  the  systea  is  governed  by  (or  X^). 
If  all  of  the  Polncar*  exponents,  have  negative  real  parts, 
or,  if  working  with  X^,  |X^|<1,  then  the  systea  is  stable. 

Usually  Ploguet  analysis  is  used  only  to  deteraine  stability, 
and  therefore  stops  at  this  point.  Por  this  study,  however,  the 
coaplete  solution  is  required,  so  P(r)  over  one  period  aust  be 
deteralned.  Substituting  equation  (2.2.11)  into  equation  (2.2. IS) 
produces 

|^(T)-A(T)P(T)-P(T)J  (2.2.17) 

Once  P(t)  is  deteralned  nuasrically,  since  it  is  periodic,  it  can 
be  expressed  in  teras  of  its  Pourier  coefficients. 

Since  both  aatrices  P  and  J  aay  be  complex,  the  above  results 
aay  be  inconvenient.  However,  they  can  be  arranged  so  that  both 
are  real.  P  should  be  arranged  in  coluan  vectors,  f^,  if  the 
eigenvalues  are  real.  If  they  are  coaplex,  the  coluan  will  be 
f.  f,.  The  aatrix  J  will  no  longer  be  in  Jordon  noraal 

vr#oii  vvfnck0 

fora,  instead  it  will  consist  of  diagonal  real  entries  <0^  if  the 

eigenvalues  are  real  or  for  a  pair  of  complex  <0^  J  will  have 

diagonal  blocks  of  the  fora; 

Held). ) 

-lm(<i)^) 


la(<i). 

\ 

Re(<i>, 


(2.2.18) 
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The  Mtrix  exp(JT)  will  contain  diagonal  entries  explu^r)  for  real 
entries/  or  diagonal  blocks; 


[co8la(w.  )t  slnln((i>.  )t 

-8lnlB((i>^)T  cosIn(a>^)T 


{2.2.19a) 


for  blocks  of  J  of  the  fora  of  Bq.  (2.2.18).  With  these  changes, 
the  previous  formulae  stay  unaltered.  Since  this  Is  a  periodic 
systea,  after  one  period  exp(<o^T)  Is  slaply  K,  therefore  the 
matrix  exp(JT)  Is  easily  constructed  as 

[Re(X.)  Ia(X.)' 

-l.(V)  Re(K;)J  <2.2.19b) 

It  Is  often  necessary,  or  at  least  aore  convenient,  to  use 
the  Inverse  aatrlx  P'^lr).  The  state  transition  matrix  Is  never 
singular  (ref  1),  therefore  F(t)  is  always  invertable.  One  method 
of  obtaining  P~^(t)  is  to  invert  F(t)  at  a  given  value  of  t. 
However,  this  method  can  be  expensive  in  both  tine  and  round  off 
errors.  Another  method  is  to  differentiate  the  Identity  FF~*si, 

^(0'“*  f  arCO-* 

‘j* 


and  the  value  for  F  substituted  from  Eg.  (2.2.17), 

|:j-|f'*J-  -F"‘(A(t)F(t)-F(t)J|  F"* 

~|f"*J«  -F~*A(t)+JF"*  (2.2.21) 

Therefore  F  (r)  can  be  constructed  by  numerical  integration  and 
transformed  into  its  Fourier  series  coefficients  without  having  to 
find  F(r).  To  develop  the  system  with  F~^(t)  the  initial 
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conditions,  or  F  (•)  nist  be  known.  Since  it  has  been  shown  frow 
equation  (2.2.131.  F(f)  is  the  aatrlx  of  eigenvectors  of  t(T) 
(which  Bay  be  coaplex),  a  real  aatrix  r(t)  can  be  constructed  in 
the  saae  Banner  as  was  developed  with  the  F  and  J  Batrices. 
Taking  the  inverse  of  this  constant  real  Bstrix. 

F‘*(i)«ri  f  T  ?  1  (2.2.22) 

I  r«al  imog  r«oi  imotg  | 
and  constructing  J  froB  the  Poincar4  exponents. 


■*(“.] 

1  • 

( 

1 

1  «•(“.] 

1  • 

( 

1 

f 

• 

Re  j 

w 

1 

K) 

6 

• 

-Ib| 

[“.] 

1  Re 

W 

F~*(r)  Bay  also  be  related  to  the  Bodal  Batrix  associated 
with  equation  (2.2.2),  i.e.,  the  adjoint  problea; 

^y(T)  «  -aVt)  (2.2.24) 

Letting  the  BOdal  Batrix  associated  with  -A*  be  L(t),  it  can  be 
found  that 

|pL(T)  .  -a’(t)L(t)+L(t)j’  (2.2.25) 

»  L(T)j’-  a'‘(t)L(t) 


or  taking  the  transpose, 

|5^^(t)  «  -l’^(t)A(t)+JL(t)’  (2.2.26) 

By  co^Mring  this  equation  with  (2.2.21)  it  can  be  seen  that 
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(2.2.27a) 


P"*(T)-L*(T) 


or 


with 


(2.2.27b) 


(2.2.27c) 


This  Illustrates  the  orthogonality  between  the  eigenvalues  o£  the 
adjoint  problem  to  the  eigenvalues  of  the  original  system. 
Solving  equation  (2.2.25)  Is  the  preferred  technique  for 
constructing  P~^(t). 

Introducing  a  new  modal  variable  n  such  that 

x(t)»P(t)7)(t)  (2.2.28) 
and  substituting  Into  equation  (2.2.2),  an  equation  for  n*  Is 
developed. 


which  reduces  to 


X*  »F* 

x'  «Ax  - ►  F'  r>^FT?'  »AF7» 

r?'  «F‘‘aFt7-F'*F*  » 

T)*  ■  j^F'*AF-F"*F*  jn 

n* «  Jf~*ap-f"*(  AF-FJ )  jn 

n'  ■  rF"*AF-F'*AF-F”*FjlT> 


7)*  •JV 


(2.2.3t) 


Therefore,  the  eigenvector  matrix  F(t)  is  a  periodic 
transformation  which  changes  the  time  periodic  system  of  equation 
(2.2.2)  to  a  constant  coefficient  system,  (2.2.31),  since  J  Is  a 
constant  matrix.  If  the  Polncar4  exponents  are  real  and  distinct, 
J  Is  diagonal  and  the  system  described  by  equation  (2.2.31) 
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c 


c 


c 


reduces  to  n  uncoupled  ordinary  differential  equations.  If  the 
Polncar*  exponents  are  coaplex,  J  is  pot  in  a  block  diagonal  fora 
with  the  real  parts  on  the  diagonal  and  the  iaaginary  parts  on  the 
off-diagonal.  Bach  coaplex  pair  results  In  a  pair  of  coupled 
equations. 


2.3  Modal  Control  Theory 

The  open-loop  characteristics  of  equation  (2.2.2)  are 

deterained  by  the  Poincar*  exponents.  These  exponents,  and 
therefore  the  systea  stability,  can  be  changed  by  adding  feedback 
control.  The  systea  has  the  fora  of  a  standard  control  problea 

xMt)«A(t)x(t)*B(t)u(t)  (2.3.1) 

where  the  control  aatrlx  B(t)  is  either  constant  of  periodic  with 

the  saae  period  T  as  the  fundaaental  systea  and  o(t)  Is  the 
control  vector.  Assuming  full  state  feedback; 

U(T)«G(T)X(T)  (2.3.2) 

where  G(r)  is  the  gain  matrix.  Denoting  the  closed-loop  state  as 
x^,  equation  (2.3.1)  I>ecoae8 

x'b(a^BG)x  (2.3.3) 

C  C 

Using  nodal  variables, 

X  »Pn 

C  C 

c  c  c 

equation  (2.3.3)  becomes 

F* »  ♦Pt>^*(A+B6)Pt7 
c  c  c 

F'rt  *Fr)  «AP77  ♦BGPk) 
c  c  c  c 


c 
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(2.3.4) 


-F‘V'  7) 

c  c  c  c 

n**  [F"*AF-P"*F*+P“‘B<!rjr)^ 

vhlch  Is  a  Ploquet  pioblea  if  It  is  assoaed  that,  at  worst,  G(t) 
is  periodic  with  period  T.  G(t)  mst  be  deterained  such  that  the 
closed-loop  systea  has  acceptable  properties.  Por  this  analysis  a 
aodal  control  technique  was  used  for  a  periodic  systea  in  which 
the  control  will  be  a  scalar  control.  This  scalar  control  was 
iapleaented  on  the  Poincar4  exponents  which  are  coaplex 
conjugates. 

Por  scalar  control  equation  (2.3.4)  is 

17'*  Jt)  ♦  g(T)u(T)  (2.3.5) 

where  g(T)  is  the  n  x  1  periodic  aode  controlablllty  aatrlx, 

g(T)  «  l’(t)  B(t)  (2.3.6) 

Modes  7)^^  are  controllable  if  the  corresponding  eleaent  g^  is 
non-zero.  In  the  case  considered  here  two  aodes  will  be 
controlled  at  once;  the  two  having  coi^ilex  conjugates  of  the 
Poincar*  exponents  with  positive  real  parts.  Therefore  the 

control  is  given  by 

u  *  lc(T)  T)  (2.3.7) 

where  k  *  (I,  f,  ...,  k^,  ...,  k^,  ...).  The  eleaents  i  and  J  aay 
be  aade  to  occur  consecutively  by  slaply  reordering  the  equations. 
The  resulting  closed-loop  equations  can  be  expressed  in  aatrlx 
fora  with  J  replaced  by  ifl  as 


n;. 
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«i>  k.g 

t 


•  "l*  ‘i«4  "u.’i 


k  a  <0  ^  k  fl 


«i> 


(2.3.8) 


Froa  equation  (2.3.8),  it  can  be  seen,  except  for  the  ith  and  i-fl 
tern,  this  cloaed-loop  systea  has  the  saae  Poincar*  exponents  as 
the  open-loop  systea.  These  two  coaponents  are  deterained  froa 
the  coupled  equations  for  the  two  aodes  and  Since  these 
two  equations  are  decoupled  froa  the  reaainder  of  equation 
(2.3.8),  they  can  be  separated  into  the  two  diaensional  systea 


k.9, 


(2.3.9) 


It  was  shorn  by  Calico  and  Viesel  (ref  1)  that,  when  g  and  g  have 

&  Z 

constant  teras  in  their  Fourier  series,  and  constant  gains  k^  and 
k^  are  used,  the  new  Poincar*  exponents  will  be  related  to  the  old 


exponents  by 


"i  ■  VV  •'s^se  (2.3.18) 
This  allows  the  sua  of  the  new  Polncar4  exponents  to  be  set  to  any 

desired  value.  While  this  does  not  provide  the  ability  to  set  the 
new  exponents  individually,  it  does  allow  the  aoveaent  of  the 
coupled  pair  of  exponents.  For  stability,  the  sua  of  these  two 
exponents  aust  be  negative.  This  is  a  necessary  and  not  a 
sufficient  condition  since  it  does  not  guarantee  the  real  parts  of 
both  exponents  are  negative  (just  their  sua).  Nuaerical 
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integration  of  the  new  aonodroMy  aatrlx  will  deternlne  the  actaal 
values  of  the  nev  Poincar4  exponents.  In  the  case  where  the 
original  Poincare  exponents  were  conplex  conjugates,  equation 
(2.3.10)  only  deteraines  the  sua  of  the  real  parts.  The  actual 
values  aay  be  either  real  or  have  iaaginary  parts. 


2.4  Developaent  of  Mon-Linear  Equations 

The  developaent  of  the  non-linear  equations  follows  froa 
Inserting  the  values  for  co  .  u  .  a  and  co  given  in  equations 

A  Z  &  Z 

(2.1.2S),  and  (2.1.42)  into  equations  (2.1.41a)  and  (2.1.41b).  As 
in  the  linear  case,  since  the  satellite  is  spinning  about  its  axis 
of  sysawtry,  is  a  constant  and  Is  a*l.  Therefore  b>^  can  be 

eliainated  froa  the  equations.  Setting  <a  in  equation  (2.1.25) 

equal  to  a-M  and  solving  for  9^  will  also  allow  the  eliaination  of 
6^  froa  the  equations. 

a+1-  ©^sln©^-  vcos©^cos©^  (2.4.1) 

Working  with  equation  (2.1.41a)  first, and  noting  that  derivatives 

with  respect  to  tiae  and  r  only  differ  by  a  factor  of  n  for  each 

order  derivative  yields, 

d  cos©  -  ©  ©  sin©  -  i^os©  sin©  +  v9  sin©  sin©  -  cos©  cos© 

(2.4.2) 

♦  vsin©^j^a+lj  ♦  |©^+  i>sin©^j©^  »  0 

Rearranging  and  further  siaplifying  yields 
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9  cos©  -  ©  ©  sin©  -  t^os©  sin©  ♦  v©  sin©  sin©^-  v9cob9  cos©^ 

ISASS  1S««CX«S 


♦  K|i©^*  ^in©  ♦  v>8in©^^  ^a+l-  ©^sin©^-  i>co8©^cos©^J=B 


©^-  ©j©^tan©j^-  vco8©^tanO^+  w©^sln©^tan©^-  i^^cos©^ 

♦  |©^+  vsin©^j  |^(K+l)(a+l)aec©^-  ©^tan©^-  vcos©^j  »  f 

©  -  ©,©  tan©  -  1.^08©  tan©.*  v6  sin©  tan©  -  i»©.cos© 
ilz  x  1  t  t  1  X  X  i 

♦  ©^(K*l)  (a*l)aec©^-  ©^©^tan©^^-  fc»©go8©^  (2.4.3) 

♦  i>(K*l)(a*l)8in©  sec©  -  u©  sin©  tan©  -v>*8in©  cos©  ■  I 

a  2  a  A  s  A  1 

This  equation  »ay  be  siaplified  and  solved  for  ©^. 

) 

♦iico8©^tan©^+  v*sin©^co8©^ 

This  is  the  sane  fora  of  the  equation  if  the  derivative  was  with 
respect  to  t.  Similar  substitutions  into  equation  (2.1.41b) 
yields, 


)  («+l)»ec©^ 

(2.4.4) 


©  »  2©^©,tan©  ♦  2v  ©  cos©  -  f©  tvs  In© 

Aa»2  S  Jb  alS 


©^+  vsin©^*i»©^cos©^-  k|©^cos©^-  vcos©^sin©^j  ^oiflj 

-  |©^cos©^-  vcos©^sln©j^ j  |a+l  -  ©^sln©^-  i^o8©j^cos©^j  (2.4.5) 

3Ksln©^cos©^ 

m  —  — ■  ■ » 

c* 

(notice  the  n*  term  is  missing  from  the  right  hand  side  trtien 
compared  to  equation  (2.1.41b)  since  the  derivative  is  now  with 
respect  to  r. 
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Slallar  steps  as  those  between  equations  (2.4.2)  and  (2.4.4)  will 
yield  froa  equation  (2.4.5)  an  expression  for  6^.  Derivatives 
with  respect  to  t  and  will  be  represented  by  dotted  teras. 


•  • 

e  « 

s 


JlC+lj  t>cosd^slnd^j  -  vslnd^-2wd^cose^co8*©^ 

-  f(9*-  p*co8  d  i^^lslnd  cosd 

m  m 


(2.4.6) 


The  expressions  for  v  and  v  can  be  found  froa  equations  (2.1.11) 
and  (2.1.12)  respectively.  The  non-linear  teras  are  now  expressed 
In  a  Banner  that  aay  be  prograaaed  and  nuaerlcally  integrated. 
These  two  equations  will  replace  the  last  two  in  equation  (2.1.54) 
and  the  feedbacic  can  be  added  In  a  slallar  aanner.  Again,  for 
'saall*  initial  conditions  this  non-linear  systea  should  behave 
siailarly  to  the  linear  systea. 
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CHAPTER  3 


RESULTS 


In  this  chapter  the  sequence  of  procedures  and  the  results  of 
the  example  tested  will  be  discussed.  The  example  was 
accomplished  employing  various  computer  programs  to  solve  the 
equations  developed  in  the  last  chapter.  The  programs  utilized 
were  run  on  a  Personal  Computer  and  were  written  In  Fortran  77 
Microsoft  version  4.fl.  This  allotted  IHSL  subroutines  to  be  used. 
Numerical  Integration  was  accomplished  by  the  NAMING  subroutine. 
It  was  decided  that  scalar  control  on  a  system  with  two  unstable 
modes  would  be  Implemented  and  the  computer  programs  were  used  to 
determine  what  parameters  would  be  used  in  the  feedback. 


3.1  Parameter  Space 

The  basic  system  was  expressed  in  the  linear  form 

X'=A(t)x  (2.2.2) 

where  the  elements  of  A(t)  are  expressed  In  equations  (2.1.54)  and 
(2.1.57).  These  elements  are  functions  of  K,  a  and  e.  Therefore 
the  stability  of  the  system  Is  determined  by  these  three 
parameters.  In  Kane  and  Barba  (4;4B5)  stability  is  shown  In  the 
K,  a  parameter  plane  for  various  values  of  the  eccentricity.  The 
values  of  the  Poincar*  exponents  for  various  values  of  K  and  oi 
with  eccentricity,  e«i.5,  is  shown  In  Table  3.1.  This  example  is 
the  same  as  used  by  Myers  (ref;5).  Myers  used  values  of  0.7,  1.0 
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ao 

3.0 

K»  1.0 

0.0000  */-  j0.1395 

0.0000  4/-  30.3525 

2.0 

0.0000  V-  j0.1174 

0.0000  4/-  30.4937 

1.0 

0.0000  */-  j0.1549 

0.0000  4/-  30.3766 

0.0 

V-  0.1878  *  j0.5000 

0.0000  4/-  30.3556  t.tt 

-1.0 

0.6718  -  0.6718 

0.0000  4/-  30.3728 

-2.0 

0.3076  -  0.3076 

0.0000  4/-  30.2952  tt 

-3.0 

V-  0.0953  ♦  0.5000 

0.0000  4/-  30.4673  t 

3.0 

K«  0.8 

0.0000  W-  j0.1489 

0.0000  4/-  30.4907 

2.0 

0.0000  t/-  j0.2241 

0.0000  4/-  30.2820 

1.0 

0.0000  */-  ii.nti 

0.0000  4/-  30.4595 

0.0 

*/-  0.0525  4  j0.5000 

0.0000  4/-  30.2638  t,tt 

-1.0 

0.6072  -  0.6072 

0.0000  4/-  30.2066 

-2.0 

0.0000  */-  j0.1272 

0.0000  4/-  30.4209 

-3.0 

4/-  0.0854  4  j0.5000 

0.0000  4/-  30.0129  t 

3.0 

K=  0.6 
0.0000  4/-  j0.1519 

0.0000  4/-  30.3791 

2.0 

0.0000  4/-  j0.0486 

0.0000  4/-  30.3082 

1.0 

4/-  0.0885  4  30.5000 

0.0000  4/-  30.2811  t,tt 

0.0 

4/-  0.0555  4  30.5000 

0.1476  -  0.1476  t 

-1.0 

0.S067  -  0.5067 

0.0000  4/-  30.0175  tt 

-2.0 

0.0000  4/-  30.1206 

0.0000  4/-  30.4005 

-3.0 

0.0000  4/-  30.4149 

0.0000  4/-  30.4414 

3.0 

K-  0.4 

0.0000  4/-  30.2559 

0.0000  4/-  30.1447 

2.0 

0.0000  4/-  30.1368 

0.0000  4/-  30.4276 

1.0 

0.0000  4/-  30.2980 

0.0000  4/-  30.0349 

0.0 

0.1249  4/-  30.2959 

-0.1249  4/-  30.2959 

-1.0 

0.2943  -  0.2943 

0.0000  4/-  30.2163 

-2.0 

0.0000  4/-  30.4390 

0.0000  4/-  30.2173 

-3.0 

0.0000  4/-  30.2863 

0.0000  4/-  30.0651 

3.0 

K=>  0.2 

0.0000  4/-  30.1215 

0.0000  4/-  30.1099 

2.0 

0.0000  4/-  30.3073 

0.0000  4/-  30.1282 

1.0 

4/-  0.0742  4  30.5000 

0.0000  4/-  30.1762  t 

0.0 

4/-  0.0710  4  30.5000 

0.0000  4/-  30.2547  t 

-1.0 

0.1860  4/-  30.4857 

-0.1860  4/-  30.4857 

-2.0 

0.0000  4/-  30.2793 

0.0000  4/-  30.4692 

-3.0 

0.0000  4/-  30.1724 

0.0000  4/-  30.4328 

Table  3.1  (ref  6)  Poincare  Exponents  for  K  vs.  a,  with  es0.5 

T  Apporvnl  prror  from  coleu  Voting  complox  logs  corroctod 
Moth  error  corrected 


3-2 


K*-B.2 


Ks-0.4 


Ks-0.6 


K»-0.8 


K»-1.0 


as 

3.0 

0.0000 

Ks  0.0 

V-  j0.0000 

0.0000 

4/-  30.0000 

2.0 

0.0000 

^/-  j0.0000 

0.0000 

4/-  30.0000 

1.0 

0.0000 

4-/-  j0.0000 

0.0000 

4/-  30.0000 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

-1.0 

0.0000 

V-  j0.0003 

0.0000 

4/-  30.0003 

tt 

-2.0 

0.0000 

4/-  j0.000« 

0.0000 

4/-  30.0000 

tt 

-3.0 

0.0000 

4/-  j0.0000 

0.0000 

4/-  30.0000 

3.0 

0.0384 

+/-  j0.0932 

-0.0384 

4/-  30.0932 

2.0 

0.0983 

+/-  J0.1756 

-0.0983 

4/-  30.1756 

1.0 

0.1829 

♦/-  j0.3833 

-0.1829 

4/-  30.3833 

0.0 

V-  0.5257 

4  j0.5000  4/ 

-0.2162 

4  30.5000 

t 

-1.0 

0.6017 

4/-  j0.0094 

-0.6017 

4/-  30.0094 

-2.0 

0.4319 

4/-  30.3953 

-0.4319 

4/-  30.3953 

-3.0 

0.2297 

4/-  30.2971 

-0.2297 

4/-  30.2971 

3.0 

0.4215 

-  0.4215 

0.0000 

4/-  30.3602 

2.0 

V-  0.4385 

4  30.5000  4/-0.3361 

4/-  30.5000 

t 

1.0 

+/-  0,6614 

4  30.5000  4/ 

-0.3876 

4/-  30.5000 

t 

0.0 

V-  0.7834 

4  30.5000  4/ 

-0.6559 

4/-  30.5000 

t 

-1.0 

0.8141 

4/-j0.0339 

-0.8141 

4/-  30.0339 

-2.0 

0.7568 

4/-j0.2810 

-0.7568 

4/-  30.2810 

-3.0 

0.4537 

4  30.4128 

-0.4537 

4/-  30.4128 

3.0 

♦/-  0.3735 

4  30.5000 

0.5384 

-  0.5384 

t 

2.0 

♦/-  0.8482 

4  30.5000 

0.5259 

-  0.5259 

t 

1.0 

f/-  1.0741 

4  30.5000  4/ 

-0.5117 

4  30.5000 

t 

0.0 

0.9275 

4/-  30.3540 

-0.9275 

4/-  30.3540 

-1.0 

0.9702 

4/-  30.0681 

-0.9702 

4/-  30.0681 

-2.0 

0.9559 

4/-  30.1550 

-0.9559 

4/-  30.1550 

-3.0 

0.8834 

4/-  30.3668 

-0.8834 

4/-  30.3668 

3.0 

♦/-  1.2424 

4  30.5000  4/ 

-0.6066 

4  30.5000 

t 

2.0 

+/-  1.2570 

4  30.5000  4/ 

-0.7543 

4  30.5000 

t 

1.0 

+/-  1.1584 

4  30.5000  4/ 

-0.9504 

4  30.5000 

t 

0.0 

1.0840 

4/-  30.2540 

-1.0840 

4/-  30.2540 

-1.0 

1.0985 

4/-  30.1107 

-1.0985 

4/-  30.1107 

-2.0 

1.1002 

4/-  30.0117 

-1.1002 

4/-  30.0117 

-3.0 

1.0896 

4/-  30.1272 

-1.0896 

4/-  30.1272 

3.0 

1.2095 

4/-  30.1618 

-1.2095 

4/-  30.1618 

2.0 

1.2095 

4/-  30.1618 

-1.2095 

4/-  30.1618 

1.0 

1.2095 

4/-  30.1618 

-1.2095 

4/-  30.1618 

0.0 

1.2095 

4/-  30.1618 

-1.2095 

4/-  30.1618 

-1.0 

1.2095 

4/-  30.1618 

-1.2095 

4/-  30.1618 

-2.0 

1.2095 

4/-  30.1618 

-1.2095 

4/-  30.1618 

-3.0 

1.2095 

4/-  30.1618 

-1.2095 

4/-  30.1618 

Table  3.1  (continued) 
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and  f.5  for  K,  a  and  e  respectively.  These  values  resulted  in  two 
conplex  conjugate  pairs  with  two  stable  and  two  unstable  roots. 

Modal  variables  are  used  to  decouple  the  four  equations  into 
two  coupled  pairs.  Phase  plots  can  then  be  used  to  view  the 
stability  and  the  effect  of  control  on  the  system.  Since  the 
modal  variables  are  four  in  number,  graphical  presentation  is 
difficult.  Therefore  phase  portraits  vill  plot  n  versus  n  and 

S  Z 

7)  versus  t)  .  In  the  non-linear  system,  since  modes  3  and  4  are 
stable,  emphasis  vill  be  place  on  modes  1  and  2  with  initial 
conditions  of  7)^  and  set  to  zero.  Modes  3  and  4  are  monitored 
to  insure  they  do  not  become  unstable,  but  the  motion  resulting 
from  the  control  to  stabilize  modes  1  and  2  are  the  primary 
concern.  Once  the  feedback  is  iBq>lemented,  the  original  modal 
variables  vill  no  longer  decouple  the  tvo  sets  of  modes.  Hovever, 
useful  Information  can  be  obtained  from  this  set  of  original  modal 
variables.  In  order  for  the  nev  equations  to  be  decoupled,  a  nev 
J  matrix  needs  to  be  developed  and  found  numerically  for  the  nev 
controlled  system.  This  vill  be  left  to  later  Investigations. 

3.2  Uncontrolled  Linearized  System 

The  parameters  used  in  this  thesis  are  K^S.?,  o=1.0  and 

e=l.5.  This  represents  the  most  common  type  of  instability,  but 
the  methods  used  to  control  this  case  are  applicable  to  other 
cases.  With  no  control  added,  the  characteristic  exponents  are; 

«  f.Sm  +/-  ji.3667  unstable 

•  a  * 

u  -9.»m  +/-  ji.3667  stable 
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Vhlle  all  the  motions  of  the  satellite  are  Important,  a  variable  ^ 
which  Is  the  angle  between  the  satellite  spin  axis  and  the  orbit 


normal,  shows  the  combined  movement  of  6  and  6  and  Is  defined  by 
^  -  arccos(cosd  *  cosB  )  (3.2.1) 

(See  Figure  (2.2))  The  motion  of  the  uncontrolled  system  Is  shown 
In  Figures  (3.1)  through  (3.11).  Figures  (3.1)  through  (3.4)  show 
the  motion  of  the  state  variables.  This  motion  Is  started  by  an 
Initial  displacement,  l.e.,  the  state  variable's  Initial  values. 
In  radians  are 


x(0) 


8.S2 

0.02 

0.00 

0.00 


(3.2.2) 


Figure  (3.5)  reveals  the  combined  motion  of  9  and  9  .  Figures 

X  z 

(3.6)  through  (3.9)  display  the  modal  response.  The  Initial 
values  for  the  modal  variables  axe  defined  by  the  use  of  Inverting 
equation  (2.2.28)  hence; 

T)(0)=F(0)’‘x(0) 

or  (3.2.3) 

T)(0)=L(0)’x(0) 

Figures  (3.10)  and  (3.11)  are  the  phase  portraits  of  the 
linearized  uncontrolled  responses.  It  can  be  seen  from  these 
modal  plots  that  the  first  mode  pair  spiral  outward  and  therefore 
Is  unstable  while  tlie  second  pair  spirals  Inward  Indicating  a 
stable  mode  pair. 


ncontrolled  Linearized  Xy  response 
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TIME  (periods) 

Uncontrolled  Linearized  X4  response 


Figure  3.6  —  Uncontrolled  Linearized  First  Mode  response 


Uncontrolled  Linearized  Second  Mode  response 


igure  3.6  —  Uncontrolled  Lir  lorized  Third  Mode  response 


Uncontrolled  Linearized  Fourth  Mode  response 


Fiqure  3.11  -  Phase  Portrait  for  Second  Pair  of  Uncontrolled  Linearized  Modes 


3.3  Scalar  Control  of  Linearized  Systen 

FroB  Floguet  theory  in  chapter  2,  a  solution  to  the  linear 

system  of  equation  (2.2.2)  is  given  by  equation  (2.2.11); 

»(T)«F(T)e'*V‘(e) 

with  J  constructed  as 

'  S.fl77  -0.3667  0.0000 

0.3667  -0.0177  0.0000 

”  0.0000  0.0000  -0.0177 

0.0000  0.0000  0.3667 

The  values  for  the  elements  of  the  periodic  F  matrix  could  have 
been  found  by  numerical  integration  and  a  computer  program  used  to 
express  these  elements  in  terms  of  their  Fourier  coefficients. 
Instead,  the  adjoint  problem  was  solved,  hence  the  values  of 
were  needed.  This  matrix  was  expressed  as  the  L(t)  matrix 
and  the  Fourier  coefficients  of  the  L  matrix  were  found.  A 
typical  element  of  L  is 

1^^=  0.9242  -  0.6855  cos(t)  -  0.4349  cos(2t)  -  0.0124  cos(3t) 

-i-  0.0011  cos(4t)  f  0.0271  cos(5t)  *  0.0354  cos(6t) 

^  0.0340  cos( /I ;  f  . . . 

+  0.6927  sin(T)  -  0.1637  sin(2T)  +  0.1957  sin(3T) 

-  0.1269  sln(4T)  -  0.0503  sin(5T)  -  0.0160  sin(6T) 

-  0.0005  sin(7T)  f  ... 

where  the  values  of  the  higher  frequency  coefficients  become 
increasingly  smaller  in  magnitude.  A  list  of  the  Fourier 

coefficients  of  selected  elements  of  the  L  matrix  is  included  in 
Appendix  A.  The  linearized  system  has  two  unstable  modes  and  the 
control  Implemented  is  angular  rate  feedbacic.  In  this  case  the 


0.0000 

0.0000 

-0.3667 

-0.0177 


(2.2.11) 


(3.3.1) 
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) 


form  of  egaation  (2.3.8)  will  be  used  vlth  the  vectors  k  and  B  as 

k  «  Jk^k^i  fj  (3.3.2) 


and 


(3.3.3) 


Prom  equation  (2.3.€)  this  produces  a 

n 


♦ 

1 

81 

41 

+ 

1 

82 

42 

♦ 

1 

88 

48 

.  + 

1 

84 

44 

9  vector  of  the  form 


(3.3.4) 


Using  the  trace  method  described  by 
gives  equation  (2.3.11)  in  the  form 

<ii'*  (O'  •  <il  *  0)  ■t'  k  il  4l  I 

IS  i  2  1  Si  4iJ 

Where  the  values  of  1  will  be  the  constant  terms  of  the  Fourier 
coefficients.  Therefore  the  first  two  elements  of  g  are  chosen 
to  be 


Calico  and  Viesel  (1;674) 

*  ‘.(v »«) 


0.1967  (3.3.b) 

g,  «  -i.8486  (3.3.7) 

The  list  of  Fourier  coefficients  for  the  first  two  elements  of  the 
g^  vector  is  also  included  in  ^pendlx  A.  As  stated  in  chapter  2, 
while  the  values  of  the  new  characteristic  exponents  can  not  be 
individually  chosen,  their  sum  can. 

<o^*(0'  »  2(0.0177)  ♦  k^(B.1967)  ♦  k^(-0.8486)  (3.3.8) 

The  values  k^  and  k^  are  picked  to  be  '-0.4  and  0.4  respectively. 


This  makes  the  new  sum 


C»'*  o>'  S  f.fJSI  -  i. 4(1. 1453)  «  -f.asi?  (3.3.9) 

Integrating  the  new  ff  matrix  over  a  period  of  2n  and  using  Ploquet 
theory  yields  the  new  characteristic  exponents, 

ta'  >  -«.f454  4  ji.SBti  and  «  -f.3373  4  jf.Sfff 
The  imaginary  parts  of  the  new  exponents  exist  because  with 

control  the  period  is  now  twice  the  old  period.  Vhen  the  new  ft  is 
integrated  over  4n,  the  values  of  the  characteristic  exponents  are 
the  real  values  with  no  Imaginary  part.  The  resulting  controlled 
motion  is  shown  in  Figures  (3-12)  through  (3-14).  Figure  (3-12) 
show  the  combined  0  and  0  response  and  Figures  (3-13)  and 

*  Z 

(3-14)  are  phase  portraits  of  the  controlled  linearized  modes 
using  the  uncontrolled  modal  variables.  Figure  (3-13)  shows  how 
fast  the  first  mode  pair  now  converges  to  a  stable  condition  while 
it  can  be  seen  in  Figure  (3-14)  the  second  pair  is  effected  only 
slightly  and  remains  stable.  For  these  plots  to  be  as  smooth  as 

Figures  (3-10)  and  (3-11)  new  modal  variables  would  have  to  be 

defined  for  the  control  system. 
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3.4  Non-Linear Ized  Systea 

With  'saall*  dlstazbances  the  uncontrolled  non-llnear  systea 
behaves,  as  predicted  by  theory,  like  the  linear  systea. 
Coaparlng  Figures  (3.15)  through  (3.21)  to  the  linearized  systea 
figures,  shows  how  close  the  uncontrolled  response  follows  the 
linearized  systea.  The  Irregular  aotlon  toward  the  end  of  the 
trace  on  figure  (3.21)  Is  due  to  the  coupling  of  the  pair 

since  they  are  Increasing  In  magnitude  as  time  gets  large.  The 

controlled  non- 11 near Ized  systea  duplicates  the  controlled 
linearized  response  since  the  coupled  effect  of  the  pair 

does  not  appear  since  the  controlled  pair  Is  now  stable.  This  Is 

Illustrated  In  Figures  (3.22)  through  (3.24). 

The  phase  portraits  of  the  controlled  non-llnear  systea  show 
the  effect  of  control  on  the  systea.  For  this  study.  Initial 
disturbances  were  selected  for  the  aodal  variables  7)^  and  but 
the  Initial  values  of  and  were  selected  to  be  zero.  The 
linearized  system  has  uncoupled  aodal  responses,  however  as 
nonlinear  effects  becoae  apparent  the  two  separate  phase  planes 
»  ,  7)  and  7?  ,  77  become  coupled.  This  is  the  reason  why,  as  the 

*  S  8  4 

disturbances  get  large,  initial  points  seeaingly  close  together 
display  very  different  responses. 

With  the  aagnltude  of  the  disturbances  for  77  and  77  kept 

»  S 

below  f.25,  the  response  Is  very  slallar  to  a  linearized  system. 
This  Is  shown  in  Figures  (3.25)  through  (3.28).  As  the 
disturbances  grow,  there  Is  an  increasing  tendency  to  oscillate 
around  the  zero  disturbance  instead  of  returning  to  complete 
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Uncontrolled  Non— Linearized  X4  response 
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ase  Portrait  of  Controlled  Non— Linearized  tji  vs.  t)2 


magnitude  of  0.25  and  same  sign 


stability.  This  Indicates  a  region  of  chaos.  This  tendency 
appears  to  Increase  sooner  in  the  second  and  fourth  quadrants. 
Figures  (3.29)  through  (3.31)  illustrate  this  tendency.  It  can  be 
seen  in  all  the  Phase  Portraits  that  the  plots  are  sysaetrical 
about  the  origin  even  in  the  non-linear  system.  With  even  larger 
initial  displacements  the  stable  areas  become  narrover  even  in  the 
first  and  third  quadrant.  The  points  shown  in  Figure  (3.32) 
without  curves  are  initial  values  which  do  not  converge  on  the 
origin  but  oscillate.  The  traces  of  the  curves  that  do  converge 
in  these  figures  appear  to  enter  the  stable  area  along  the 
vertical  axis. 

If  the  new  modal  variables  \nze  found  for  the  controlled 
system  these  plots  should  become  smoother  with  the  lines  not 
crossing  in  the  linear  region  since  these  new  modal  variables 
would  again  decouple  the  the  two  mode  pairs.  Using  the 
uncontrolled  modal  variables  still  allows  conclusions  to  be  made 
with  respect  to  the  size  of  'small'  disturbances.  It  is  evident 
from  the  phase  portraits  that  a  modal  displacement  magnitude  less 
than  0.25  radians  and  radians/r  will  result  in  a  stable  system 
with  reducing  magnitudes.  Larger  displacements  will  result  in  a 
reducing  magnitude  or  chaos.  Further  investigation  is  warranted 
in  this  area  to  further  define  the  value  for  'small*.  The  use  of 
controlled  modal  variables  will  produce  more  Information. 
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Chapter  4 

Conclusions  and  Reconanendations 

The  results  of  the  previous  chapter  leads  to  sone  important 
conclusions  on  the  concept  of  small  disturbances  of  a  satellite 
system  and  the  resulting  motion  of  that  satellite.  It  was  shown 
scalar  control  works  reasonability  well  on  the  system  with  two 
stable  and  two  unstable  modes.  Left  to  be  determined  is  a  method 
in  selecting  the  optimum  gain  of  the  feedback  other  than  'hit  and 
miss'.  Large  gains  do  not  necessarily  produce  the  most  stable 
system.  While  large  gains  may  produce  one  root  which  is  more 
stable,  it  is  likely  to  push  another  root  in  the  opposite 
direction  toward  instability.  The  modal  variables  used  in  this 
thesis  is  for  the  uncontrolled  system  and  therefore  do  not 
decouple  the  controlled  modes.  However,  the  uncontrolled  modal 
variables  can  provide  information  on  the  stability  movement  of  the 
controlled  modes,  but  further  integration  to  obtain  the  new 
controlled  modal  variables  would  provide  better  information. 

From  the  uncontrolled  modal  variables  it  is  apparent  that 
there  are  several  regions  of  different  stability.  Under  a  value 
of  0.25  the  non-linear  equations  behave  linear.  Outside  this 
region,  there  exists  an  area  where  different  initial  conditions 
may  or  may  not  return  the  system  to  its  original  equilibrium 
point.  Starting  points,  seemingly  close  together  on  a  two 
dimensional  plot,  result  in  very  different  reactions.  Some  return 
to  the  center  while  others  seem  to  oscillate  in  undamped  but 
bounded  movement  about  the  center.  This  can  be  attributed  to  the 
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two  diaensional  representation  of  a  four  dinensional  problem  and 
to  the  use  of  uncontrolled  modal  variables  In  the  controlled 
system. 

Further  study  Is  warranted  In  four  areas.  First,  the  new 
modal  variables  should  be  found  to  further  this  Investigation. 
Then  this  outer  area  of  confusion  can  be  studied.  Also  a  method 
needs  to  be  found  for  systematically  determining  the  optimum  gain 
for  the  feedback.  Finally,  study  of  control  other  than  scalar 
should  be  Investigated  using  the  non-llnear  equations.  Work  on 
vector  control  with  the  controlled  modal  variables  seems  to  be  the 
next  logical  step. 
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Appendix  A:  Foariex  Coefficients 
Adjoint  Solution  Hatrixx  (first  10  terns) 


n 

a 

n 

b 

n 

II 

.92423111685947100^00 

H 

-.68545295090251760^00 

.6926555705680942400 

B 

-.43488643796333080400 

-.1636572707973497400 

B 

-.12435058525208430-01 

.1957220216151334400 

B 

.11136936380036850-02 

-.1268693977398193400 

5 

.27146253411323960-01 

-.5032345220286653-01 

6 

.35366289204882160-01 

-.1604309578632710-01 

B 

.33984111360418210-01 

-.4595563214445126-03 

B 

.28771080265566940-01 

.5652832193369086-02 

9 

.22727932659096520-01 

.7218445776603422-02 

10 

.17191192946061540-01 

.6801233935661503-02 

n 

a 

n 

b 

n 

0 

.1225952800437708B401 

1 

-.22124209ie689654B401 

.58202064134e0934B400 

2 

-.7401668126789335E400 

.4269529560564276E400 

3 

-.5228195797654754E400 

.6112e75661240712E-01 

4 

-.6759642857946536B-02 

.1662875661240712B-S3 

5 

.2060223B91165860E-01 

-.2817629049667971B-01 

6 

.3186080763242387E-01 

-.3603346602853508B-01 

7 

.3169692364656334B-01 

-.3422604876598152B-01 

8 

.2713412076597691K-01 

-.2876258140367391B-01 

9 

.2152152407587097E-01 

-.2260715215954215E-01 

.1630200312669608E-01 

-.1703811298128566B-01 
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n  a 

n 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


.1967460398576639B400 

.5173699814306124B400 

.7676230805236048B-02 

.9436071392305762B^00 

-.2651106267160906B400 

-.2170659379275563B-01 

-.2406769170357917B-02 

.1101699016076314B-02 

.1505842531590358B-02 

.1237808414617844B-02 

.8913929114430876B-03 

.6068758572467093B-03 

.4013494912402667B-03 

.2610136202602885B-03 

.1679924698986054B-03 

.10739458883868b6E-03 

.6834469788277933B-04 

.4335866054233590B-04 

.2744787174179565B-04 

.1734952804595234B-04 

.1095500361219531B-04 

.6912323480100846B-05 

.4359319069900982E-05 

.2748242086847370B-05 

.1732024486001318B-05 

.1091170288716139E-05 

.6870500555549305E-06 

.4321902615741710E-06 

.2714311273037677B-06 

.1700032999699712B-06 

.1059947298799529E-06 


(first  30  terns) 
b 

n 


.1611886143123140B400 

.6897673595227186B400 

.6599075565223430E1^00 

-.1288616102659364Bf00 

.5790067376471739E-01 

.4576824322517970E-01 

.3051217699305379E-01 

.1953038243976544E-01 

.1231450431329914B-01 

.7714018512527845B-02 

.4817728885472591E-02 

.3004941590834503E-02 

.1873444969883763E-02 

.1168053395432101B-02 

.7284763322974062E-03 

.4545299879918495B-03 

.2837514607110684E-03 

.1772389307880157E-03 

.1107726452879897E-03 

.6927229730155386E-04 

.4334497409023176E-04 

.2713735810857480E-04 

.1699981558207768B-04 

.1065537202337215B-04 

.6682569071913408E-05 

.4193482545194881E-05 

.2633131743897344E-05 

.1654434650785011E-05 

.1040206602273654B-05 

.6544726496219521E-06 
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